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ABSTRACT 

Over the last few years increasing consideration has been given to the study of Laser 
Guide Stars (LGS) for the measurement of the disturbance introduced by the at- 
mosphere in optical and near-infrared astronomical observations from the ground. A 
possible method for the generation of a LGS is the excitation of the Sodium layer in 
the upper atmosphere at approximately 90 km of altitude. Since the Sodium layer is 
approximately 10 km thick, the artificial reference source looks elongated, especially 
when observed from the edge of a large aperture. The spot elongation strongly limits 
the performance of the most common wavefront sensors. The centroiding accuracy 
in a Shack-Hartmann wavefront sensor, for instance, decreases proportionally to the 
elongation (in a photon noise dominated regime). To compensate for this effect a 
straightforward solution is to increase the laser power, i.e. to increase the number of 
detected photons per subaperture. The scope of the work presented in this paper is 
twofold: an analysis of the performance of the Weighted Center of Gravity algorithm 
for centroiding with elongated spots and the determination of the required number of 
photons to achieve a certain average wavefront error over the telescope aperture. 

Key words: instrumentation: adaptive optics - atmospheric effects - telescope 
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1 INTRODUCTION 

Adaptive Optics (AO) systems may provide diffraction- 
limited imaging with ground-based telescopes, measuring 
the wavefront aberrations due to the atmosphere by means 
of a suitably bright reference star and compensating for 
these aberrations by means of a deformable mirror. In clas- 
sical AO systems, based on a single natural reference star, 
such a star has to be located close to the target object, i.e. 
within the isoplanatic patch, so that the probability of find- 
ing it is quite low. This probability defines the fraction of 
the sky where an efiicient correction can be achieved (sky 
cover age), and for near infrared wavelengths it is few per- 
cent (|Beckerj|l993 ). In this context the po tential power of 
Laser Guide Star (LGS) AO, proposed bv IFov fc Labevrig 
l| 19851 ). is clear. The basic idea is to generate a bright 
guide star in the upper atmosphere by means of a laser. 
Two alternative physical phenomena may be exploited for 
such a purpose: Rayleigh scat tering at low-medi um alti- 
tude, up to ~ 10-20 km (jFugate et all Il99ll ). or ex- 
citation of the atoms in the Sodium layer located at ap- 



proximately H ~ 90 km, characterized by a thickness of 
/^H ^ 10 km ijThompson fc Gardner! 1 19871 : [Gardner et~al] 
1 19901 : iHapper et al.lll994l ). Although low-altitude Rayleigh 
scattering can provide much more backscattered light, the 
higher altit ude of the mesospheric Sodium layer reduces the 
cone effect (|Fov fc Labevrielll985l ). i.e. maximizes the atmo- 
spheric volume sampled by the artificial source. 

The use of LGS in principle solves the problem of 
the lack of bright reference sources, but other difficul- 
ties arise. Apart from the above mentioned cone effect, 
that might be compens ated using multiple reference sources 
(|Ragazzoni et al]| 19991 ). the problem analyzed in this paper 
is related to the perspective elongation of the LGS image, 
due to the vertical extension of the artificial source. Refer- 
ring, for instance, to a Shack-Hartmann (SH) WaveFront 
Sensor (WFS), the subapertures located far from the laser 
launcher see an elongated LGS image, pointing towards the 
laser launcher. The angular elongation is given by (see Fig- 
ure [TJ 
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Figure 1. Elongation of the LGS image at a distance L (in this 
case at the edge of the telescope aperture) from the launch site. 
l\H represents the thickness of the sodium layer, H is the mean 
altitude of the layer and O^^ia is the angular elongation. 
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Since a and (5 axe given by 
Equation ([T}, in the small angles approximation, becomes 



m + HAH 



(2) 



where AH is the thickness of the sodium layer, H is the 
mean altitude of the sodium layer and L is the distance be- 
tween the pupil subaperture and the LGS launcher (Figure 
[T]). This elongation can dramatically affect the capability to 
measure the wavefront, especially in the case of large diam- 
eter telescopes. 

On a 30-40 m class teles cope, hke the TMT 
llNelson et al.ll2006l ) and the E-ELT l|Gilmozzi fc Spvromilicl 
l2008l ). the maximum elongation 9e.io varies between few arc- 
sec and 10 arcsec, depending on the actual telescope diame- 
ter, on the Sodium layer properties and on the laser launcher 
position. 

Through numerical simulations, verified by means of an- 
alytical formulas, we compute the measurement error on the 
Optical Path Difference (OPD) per subaperture, consider- 
ing t he elongation pattern on a SH WFS for the 42 m E- 
ELT l|Gilmozzi fc Spvromiliol bOOSl ). in the case of a laser 
guide star launched either from behind the secondary mir- 
ror (i.e. from the center of the pupil) or from the edge of 
the primary mirror. In Section [5] we describe the simula- 
tions and the input parameters. The performance analysis 
of a specific algorithm (Weighted Center of Gravity) is the 
subject of Section [31 In Sections [3] and [S] we propagate the 
measurement error on the OPD per subaperture to RMS 
wavefront error over the whole pupil, hereafter referred to 
as WFE, assuming a modal least square wavefront recon- 
struction on a Karhunen-Loeve (K-L) basis. The resulting 
WFE is evaluated as a function of number of photons and 
Read-Out-Noise (RON). An estimation of the number of de- 
tected photons needed to achieve a pre-fixed WFE per LGS 
is given in Section (5] 



Figure 3. Simulated LGS images with different elongations, as 
seen from subapertures at different distances from the laser pro- 
jector, considering a Gaussian profile (top panel) and a bi-modal 
density profile (bottom panel). 



2 METHOD 

2.1 Sodium density profile 

In order to understand how the image of a LGS looks like as 
seen from each subaperture of a Shack-Hartmann wavefront 
sensor, it is important to choose a model of the Sodium 
density profile. The density of the Sodium layer can change 
by a factor of two or more depending on the telescope site 
and the time of the y ear. Its mean altitude can change on 
a time scale of hours (jO'SuUivan et al1l2000l ). We consider 
two different Sodium profile models (Figure [H), with a shape 
similar to those presented in lDrummond et al.l (|2004h : 

• Single Gaussian Sodium density profile, peak at // = 90 
km and FWHM AH = 10 km; 

• Bi-modal Sodium density profile given by the sum of 
two Gaussians with the following characteristics: peak at 
// = 84 km and FWHM AH = 8.24 km, peak a.t H = 94.5 
km and FWHM AH = 2.35 km. 



2.2 Modeling elongated spots 

The described Sodium profiles have been used to simulate 
the spots seen by a SH WFS. Some examples of simulated 
spots are shown in Figure [3] 

We consider two LGS launching schemes: from behind 
the secondary mirror (central projection) and from the edge 
of the primary mirror (lateral projection) . Central projection 
gives a perspective elongation pattern with radial symmetry, 
while lateral projection leads to a non-symmetric pattern 
with maximum elongation two times larger than in the for- 
mer option. The peak density of the Sodium layer is imaged 
at the center of the Field of View (FoV) of each subaperture. 
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The LGS spot imaged by each subaperture is modeled 
as follows: 

• the starting point is a line source of length calculated 
accordingly to Equation [2] and with intensity distribution 
given by the adopted Sodium density profile; 

• the line source is convolved by a seeing-like disk of an- 
gular size 9s — 1.5 arcsec, taking into account the round-trip 
propagation of the beam through the turbulent atmosphere 
(seeing = 1 arcsec). This angular size corresponds to the 
width in the non-elongated axis. 

The subaperture diffraction {Od 0.24 arcsec) is not consid- 
ered in the spot model, as it is negligible with respect to the 
seeing effect, that always dominates over diffraction since 
the AO correction at the laser wavelength is quite poor. 

We generate sequences of instantaneous spots: each spot 
is characterized by a random jitter and is contaminated by 
photon noise and read-out noise (RON). The random offset 
is extracted from a Gaussian distribution, whose standard 
deviation may be tuned to simulate open or closed loop con- 
ditions. The photon noise is added to each pixel accordingly 
to a Poisson distribution with a mean value equal to the 
pixel intensity; the RON is extracted from a Gaussian dis- 
tribution of zero mean. 

Different values of focal plane sampling and subaper- 
ture FoV may be adopted in the generation of the sequence, 
in order to evaluate their impact on the WFS performance. 
Depending on the chosen subaperture FoV, the wings of the 
most elongated spots, especially for the the subapertures 
at the edge of the pupil, may fall outside the FoV of the 
subapertures themselves. In a real system, this would trans- 
late into a light contamination of adjacent subapertures, a 
drawback that may be overcome using simple light baffling 
methods, for instance a diaphragm on the LGS image in 
an intermediate focal plane before the WFS. This baffling 
is always implicitly assumed here, so that the limited sub- 
aperture FoV translates into a truncation of the spot wings, 
but not into a contamination of adjacent subapertures. The 
number of photons per subaperture referred to in the fol- 
lowing includes also the photons blocked by the subaperture 
baffle. 



2.3 Algorithm performance evaluation 

The sequence of instantaneous spots, typically formed by 
500 realizations, is analyzed by a spot position measure- 
ment algorithm. The algorithm result is compared to the 
known true position, thus obtaining a sequence of position 
differences, the RMS of which is the figure of merit for the 
algorithm performance evaluation. Each instantaneous spot 
is tested for detectability. For thi s purpose two checks are 
implemented l|Thomas et al.|[2006l '): the maximum intensity 
value of the spot is requested to be at least twice the RON 
and the calculated position has to be within 1 — a from the 
true position, where cr is the standard deviation of the best 
fit Gaussian to the long-exposure spot, obtained by sum- 
ming all the images in the sequence. When more than a 
certain fraction of images (in our simulations 50%, following 
Thomas et al. 2006) are rejected, we consider the centroid 
measurements to have failed and not be reliable for those 
light conditions. 

The main simulation parameters are listed in Table [1] 



Telescope parameters 


D = 42 m 




Relative obstruction 


0.3 linear 




M 1 J. J ■ TTITI 7'T TTV T 

Non-elongated axis b WHJVL 


~ 1.5 arcsec 




Single Gaussian Profile 


iU km t WrlM '& 


90 km 


Double Gaussian Prolilc 


8.24 km FWHM 


@84 km 




2.6o km t WrlM 


(Q94.5 km 


No. of subapertures across D 


84 




Subaperture FoV 


12x12 pixel 




Pixel scale 


0.75 arcsec/pixel 




No. of photons per subap. 


up to 2000 




RON 


1, 3, 5, 7, 9, 11 e 


~ / pixel 


Number of K-L modes 


5600 





Table 1. Simulation parameters. 



The values of some parameters, like the number of pixels 
per subaperture and the pixel scale, are the result of an 
optimization presented in Section 13.21 The performance of 
the algorithm has been evaluated considering subapertures 
at difference distances from the laser launcher, in order to 
have a complete map across the whole telescope pupil. 



3 THE WEIGHTED CENTER OF GRAVITY 
ALGORITHM 

3.1 Algorithm description 

The estimation of the position of the spot on the focal plane 
of each subaperture of the SLI WFS is generally performed 
by computing the Center of Gravity (CoG) of the spot. In 
this paper we analyze the performance of a slightly different 
algorith m, the Weigh t ed C enter of Gravity (WCoG), pro- 
posed bv lFusco et al.l (j2004l ). This algorithm is a maximum 
likelihood estimator of the spot position in presence of Gaus- 
sian noise. The expression of the x component of the spot 
position in a cartesian {x,y) reference frame is given by 



Cx 



(3) 



where Xi^j is the x coordinate of pixel (i,i), is the in- 
tensity of pixel (i, j) and Wij is a suitable weighting func- 
tion. The y component of the spot position is si milar. The 
theoretical weighting function (|Fusco et al.ll2004h is the in- 
stantaneous noiseless spot image itself. In the case of a small 
spot jitter in the subaperture, the weighting function can be 
assumed equal to the mean spot, i.e. to the average of the 
500 realizations in our case. 

The weighting function allows the attenuation of the 
noise effects, especially in the pixels with a low signal, but it 
also introduces a distortion on the position estimation, pro- 
portional to the distance of the actu al spot position from the 
center of the weighting function. In iThomas et all (|2006l ). a 
correction factor, called 7 factor, has been analytically calcu- 
lated in the case of a Gaussian spot. Since we are interested 
also in bi-modal and more general spot shapes and since 
a good knowledge of the instantaneous spot and weighting 
funct ion are crucial to calculate the 7 factor (jThomas et al.l 
12008"), instead of calculating a multiplicative 7 factor, we em- 
pirically derived a calibration curve to compensate for the 
biasing effects due to the weighting function. For a given 
spot template, depending on the Sodium profile, the result 
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Analyticol formulas (Nsomp — 0,3) 

Simulotion results Nph=250 



-0.6 -0.4 -0.2 -0.0 0.2 0.4 0.6 

Actual spot position (pixel) 

Figure 4. Example of calibration curve for the WCoG algorithm. 
A template is moved in a set of known positions and then the 
spot position is measured by mean of the WCoG. The estimated 
positions are plotted against the actual ones. The resulting plot is 
fitted by a polynomial function, that is used to correct the WCoG 
measurement for any position inside the subaperture FoV. 



of the WCoG applied to the template not contaminated with 
noise is plotted against the actual spot position for several 
known positions inside the subaperture FoV. The resulting 
plot is the calibration curve that allows the compensation of 
the bias. The spot template can be thought of as the average 
spot over a period of time over which the Sodium profile is 
approximately constant. The calibration curve depends on 
the spot elongation and therefore changes with the subaper- 
ture position in the pupil. For each subaperture we find a 
curve like the one shown in Figure |4] The calibration curve 
is fitted by a second order polynomial, which is used to re- 
trieve the correct spot position for any position inside the 
FoV. This method easily allows the compensation of non- 
linear effects, due for instance to the truncation of the spot 
at the subaperture FoV edge or to sampling effects like in 
a quad- cell scheme, a s discussed in a paper in preparation 
l|Schreiber et allbOOSD . 

The results of the numerical sim ulations have b e en ver - 
ified using the analytical formulas of lThomas et al.l l|2008h : 
the parameter Nsamp of that paper, defined as the number 
of pixels per FWHM of the subaperture diffraction limited 
PSF, is assigned here the value Nsamp ~ 0.3 (Figure [5]). 

A good agreement has been found between the simula- 
tions and the analytical formulas. A slight disagreement is 
found in the case of low number of photons, probably be- 
cause of the increased measurement noise due to the pres- 
ence of negative pixels. 



3.2 Algorithm optimization and performance 

The measurement error on the OPD per subaperture is 
taken as a figure of merit for the optimization of the algo- 
rithm parameters and for the performance evaluation. The 
basic parameters are the focal plane sampling and the FoV of 
the subaperture. Given a subaperture FoV, the calculation 
may be restricted to a smaller window, to optimize the per- 




12 3 4 5 6 

elongated axis FWHM [orcsec] 

Figure 5. Comparison between OPD measurement error as a 
function of the elongation value, parametric in the number of 
photons, calculated by numerical simulations (symbols) and by 
the analytical formula (continuous line). Central laser projection, 
RON = 3 e- /pixel. 



formance on particular Sodium profiles, like the bi-modal, as 
will be discussed later in this Section. For the optimization 
of the parameters, a reasonably high signal-to-noise ratio 
case has been considered, corresponding to Nph =500 and 
RON — 3 e~ /pixel. A central laser projection scheme has 
been assumed. 

The first step of the optimization concerns itself with 
the sampling. Assuming a wide FoV, wide enough not to 
introduce any relevant truncation effect, the performance 
for three cases of sampling has been evaluted: 0.75, 1 and 
1.5 arcsec/pixel. Considering the non elongated spot width 
(FWHM — 1.5 arcsec, as indicated in Table[T|, the first case 
corresponds approximately to the Nyquist sampling of the 
spot itself. The performance has been evaluated in terms of 
error on the OPD measurement in an edge subaperture, as 
described in Section 12.31 for a Gaussian Sodium profile for 
different simulation sets, corresponding to different values 
of the spot jitter, ra nging from c losed to open loop condi- 
tions. Accordingly to lTvsor] l| 19981 ). the maximum spot jitter 
of 0.5 arcsec that we have considered a seeing condition of 
FWHM ~ 1 arcsec. The results are shown in Figure[Sl where 
the three upper curves refer to the elongated direction, while 
the three lower ones refer to the non elongated axis. From 
the elongated axis, one might infer that the coarser sampling 
is the best, because it improves the signal-to-noise ratio per 
pixel. However it is clear that the non elongated direction 
may be affected by undersampling effects for the coarser 
pixel scale, which increases the measurement error for the 
larger spot jitter cases. For this reason, in the following anal- 
ysis the sampling of 0.75 arcsec/pixel has been adopted in 
the following, as it assures a constant error even in the worst 
seeing conditions. We do, however, just mention that in or- 
der to slightly relax the photon ffux requirements without 
relevant drawbacks, a coarser sampling of 1 arsec/pixel also 
looks very attractive. 

The second step of the optimization is concerned with 
the subaperture FoV. The analysis has been performed on 
the two reference Sodium profiles (single Gaussian and bi- 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 

spot jitter (arcseconds) 

Figure 6. OPD measurement error in an edge subaperture for 
different values of spot jitter and for a Gaussian Sodium pro- 
file. Three different pixel scales are compared: 0.75 arcsec/pixcl 
(solid lines), 1 arcsec/pixel (dot-dashed lines) and 1.5 arcsec/pixcl 
(dashed lines). The three upper curves refer to the elongated axis, 
the three lower to the non elongated one. 
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WCoG — Gaussian Sodium Profile 
0.8 1 " I " ' " ' ^ 




0.0 I , I , , , \ , , , I , , 1 \ , , , \ , , , I , , 1 \ , 

8 10 12 14 16 18 20 

subaperture distance from the center (meters) 

Figure 7. OPD measurement error as a function of the sub- 
aperture distance from the center of the pupil (central laser pro- 
jection) for different cases of subaperture FoV. Pixel scale 0.75 
arcsec/pixel. Gaussian Sodium profile. 



modal), assuming the sampling previously derived (0.75 arc- 
sec/pixel). The OPD measurement error has been evaluated 
as a function of the subaperture distance from the laser pro- 
jector, considering different subaperture FoV values, ranging 
from 4 X 4 to 20 X 20 pixels. The results for the Gaussian 
Sodium profile and central laser projection are shown in Fig- 
ure [T] The optimal FoV for this profile is of the order of 
12 X 12 pixels, since wider FoVs do not further reduce the 
measurement error. For lateral laser projection, the optimal 
FoV is of the order of 24 x 24 pixels. The results for the bi- 
modal Sodium profile are shown in Figure [S] In this case a 
smaller FoV is clearly advantageous: this may be explained 
considering that a FoV as small as 6 x 6 pixels, centered on 
the sharper component of the bi-modal spot, translates into 
a mitigation of the spot elongation, although at the price of 
a loss of photon fiux. This analysis shows that, for central 
laser projection, the optimal subaperture FoV, accounting 
for the different Sodium profiles considered here, is of the 
order of 12 x 12 pixels, but the calculation for the bi-modal 
case is more effective if restricted to a smaller window, ex- 
tracted from the subaperture FoV. 

The WCoG algorithm performance has been then eval- 
uated for different values of photon flux and RON with the 
parameters optimized before (0.75 arcsec/pixel sampling, 
12 x 12 pixels FoV with windowing down to 6 x 6 for bi- 
modal profile). The results are shown in Figures [9l and [10] 



4 NOISE PROPAGATION 

The wavefront on the whole pupil is reconstructed by means 
of the OPD measurements per subaperture and hence the 
noise associated to the measurements is propagated to the 
reconstructed wavefront. The aim of this section is the trans- 
lation of the OPD measurement errors calculated in Section 
[3] into RMS wavefront error (WFE). This propagation is 
performed assuming that the wavefront may be expanded 



WCoG - Bi-modal Sodium Profile 
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Figure 8. Same as[7l but for a bi-modal Sodium profile. The FoV 
is centered on the sharper component of the spot. 



as a linear combination of K-L polynomials; the considered 
number of modes, = 5600, corresponds approximately to 
the total number of subapertures on the pupil of a telescope 
with the specifications listed in Table [T] 

The OPD measurement problem and the wavefront re- 
constructio n process can be expressed in a matrix-algebra 
framework (|Roussetll2004l ). The measurement equation is 

s = Da + n (4) 

where a is a vector gathering the A'' unknown coeffi- 
cients of the K-L modes, s contains a set of 2M OPD mea- 
surements (M is the number of subapertures) and n is a 
vector of measurement errors. The matrix D is a 2M x N 
interaction matrix giving the wavefront sensor response to 



® 0000 RAS, MNRAS 000, 000-000 



6 L. Schreiber et al. 



WCoG — Gaussian Sodium Profile 



nph = [ 100, 250, 500, 750, 1000, 2000] 
ron = 3 
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WCoG - Bi-Modal Sodium Profile 



nph = [ 100, 250, 500, 750, 1000, 2000] 
ron = 3 




5 10 15 20 

subaperture distance from the center (meters) 



Figure 9. OPD measurement error as a function of the distance of the subaperture from the pupil center, for central laser projection. 
The OPD error is parametric in terms of the number of photons per subaperture, ranging from Nph =100 to Nph =2000. The 
lowest curve in each plot refers to the highest number of photons. RON = 3 e~/pixel, sampling 0.75 arcsec/pixel. Left: Gaussian 
Sodium profile, subaperture FoV = 12x12 pixels. Right: bi-modal Sodium profile, window size 6x6 pixel. 



WCoG — Gaussian Sodium Profile 



WGoG - Bi-Modal Sodium Profile 



nph = 500 

ron = [1 ,3,5,7,9,1 1 ] 
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Figure 10. Same as Figure|9] but curves are parametric in terms of RON, ranging from RON = 1 e /pixel to RON = 11 e /pixel. 
The lowest curve in each plot refers to the smallest value of RON. 



tlie K-L modes; apart from a normalization constant, its n- 
th column contains the average value of the derivative of the 
n-th K-L mode for each subaperture in one direction. The 
matrix D is assumed to be sorted alternating the x and y 
derivatives, although any other ordering is acceptable, pro- 
vided the vector s is ordered accordingly. 

We consider a standard least-square reconstructor, 
meaning that the reconstruction matrix is simply the gener- 
alized inverse of D. Hence from Equation (|4| we can write 



where B is given by 

B = (d'^d) ^D'^. (6) 

If the measurement noise n has a covariance matrix Cs, 
the covariance matrix of the K-L coefficients is BCsB'^ and 
the wavefront error variance is the sum of the variances of 
the modal coefficients: 



aest = Bs 



(5) 



al = trace (^BC^B"^^ . (7) 
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Figure 11. Elongation pattern for a laser guide star projected 
from behind the telescope secondary mirror. (£, y) is a reference 
frame defined by the spot elongation direction and obtained by 
rotating the frame {x, y) by an angle 6. 

The noise covariance matrix Cs is a block matrix of size 
2M X 2M and can be easily constructed with the following 
prescriptions. The measurement noise in different subaper- 
tures is statistically uncorrelated, so the covariance matrix 
is a block diagonal matrix. Furthermore, if the measurement 
of the spot position could be performed with respect to the 
principal axes {x,y) of the elongated spot (Figure [TT|) . then 
the measurement noise in the two axes would be uncorre- 
lated. However the spot position measurement is performed 
with respect to the axes of the {x,y) cartesian reference 
frame and the measurements errors in x and y are generally 
correlated: they may be retrieved from the measurement er- 
rors in the (x, y) reference frame applying a rotation matrix. 
The diagonal block referring to the j-th subaperture in the 
covariance matrix Ca has the following structure: 

22,2l2 2 7 2 L\ 

■ a J + Oy. ■ bj a^^ ■ aj ■ bj ~ a^^ ■ a, ■ bj \ 

2 J. 2 u 2 1.2 , 2 2 I 

(8) 
where 

aj = cos 6j 
bj = sin Qj 

and Qj is the angle between the pre-defined global reference 
frame (x^y) and the local reference frame (x,y) defined in 
each subaperture by the direction of the elongation. An ex- 
ample of the elongation pattern on the whole aperture is 
shown in Figure 1111 

The calculation of the WFE, i.e. the square root of the 
wavefront error variance defined by Equation [T] requires the 
evaluation of very large matrices with the adopted number 
of modes (N = 5600). For the design and optimization of the 
WFS, this evaluation has to be repeated for several cases of 
RON and of number of detected photons per subaperture. 
In order to speed up the process, we have applied a known 
property o f the error pr opagation in a modal reconstruction 
approach (|Hardvlll99^ '): the propagated WFE depends on 
the number of modes N following a logarithmic law 

WFE(Af) =a + 61ogiV (9) 

where a and b are constants. The calculation of the WFE 



1000 2000 3000 4000 5000 

K-L modes 

Figure 12. Error propagation coefficient in case of no central ob- 
struction as a function of the number of modes for a case without 
spot elongation (diagonal measurement noise covariance matrix). 
Triangles represent some of the results of the calculations up to 
N = 465 modes; the extrapolated curve is the best fit logarithmic 
law described in the text. 



has been thus performed following Equation[7]for a relatively 
small number of modes, up to A'' = 465, and the results have 
been fitted by the logarithmic curve described by Equation 
[9] the best fit curve has been then extrapolated to the de- 
sired number of modes (N = 5600), representing essentially 
the maximum number of independent modes measurable 
with the telescope specified in Table [T] This method has 
been verified on a system with just 20 subapertures across 
the diameter: an excellent match has been found between 
the full computation of the WFE following Equation [7] and 
the approximated result obtained by the extrapolation of 
the best fit logarithmic curve. 

The method proposed here has been validated also in 
a limiting, though relevant, situation: the case without spot 
elongation, where the OPD measurement error is the same 
in all subapertures. Here the measurement noise covariance 
matrix is diagonal, Cs = cr^I, and the propagated error is 

a\ = trace (bb'^^ al (10) 

where cr^ denotes the measurement noise variance and 
the term trace(BB''^) is the so-called error propagation co- 
efficient, usually of the order of or lower than 1. Figure 
1121 shows the behavior of the error propagation coefficient 
for a sample of mode numbers (up to TV = 465) and the 
best fit logarithmic curve (Equation [9]), extrapolated up 
to A = 5600 modes. The error propagation coefficient es- 
timated by this extrapolation is in very good agreement 
with the value found fo llowing the method proposed by 
iRigaut fc GendronI l|l992D . 



5 REQUIRED NUMBER OF DETECTED 
PHOTONS 

Applying the method explained in the previous Section, 
we calculate the WFE in a number of cases with differ- 
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Figure 13. Wavefront error vs number of detected photons per subaperture per frame for a Gaussian Sodium profile. The 
laser is projected from the center of the pupil (left panel) and from the edge (right panel). Different curves indicates different 
values of RON. The uppermost curve in each plot refers to the highest RON considered, i.e. RON = 11 e~/pixel, and, in 
decreasing order, we find the curves relative to RON = 9, 7, 5, 3 and 1 e~ /pixel. The vertical lines represent the cut due 
to the detectability test described in Section 12.31 



WCoG 



modal Sodium Profile 




1000 1500 
photons per subaperture 



Figure 14. Same as Figure [T3l for a bi-modal Sodium profile and 
considering only central laser projection. A 6 X 6 pixels window 
has been applied to the bi-modal spot, in order to select only the 
sharper component. 



ent detected photons per subaperture per frame and dif- 
ferent values of RON. Fitting the results for a given value 
of RON with a second order polynomial curve, we find 
the WFE as a continuous function of the number of de- 
tected photons. To evaluate the required number of detected 
photons per subaperture per frame, a target value for the 
WFE has to be defined. We b ased our study on the case of 
MAORY (iDiolaiti et al.ll2008l) . the Multi-Conjugate Adap- 
tive O ptics module for the European Extremely Large Tele- 
scope (IGilmozzi fc Spvromiholliooil ). In the error budget of 
MAORY, the target for the WFE in closed loop is 100 nm. 
Considering the number of deformable mirrors, the number 
of LGSs, the noise reduction due to the closed loop correc- 



tion, the target open loop WFE applicable to the present 
single channel analysis is 200 nm. 

The analysis on the required number of photons has 
been performed in three relevant cases: 

• LGS projected from behind the telescope secondary 
mirror (maximum spot elongation 9eio ~ 5 arcsec) 

• LGS projected from the edge of the primary mirror 
(maximum spot elongation Oeio ~ 10 arcsec) 

• No elong ation, a condition achievable in principle by 
range gating (iFugate et al.| 199ll ) or dynamic refocus of a 
pulsed laser (jGeorges et al.ll2008l ). 

Figure [T3] shows the results for the Gaussian profile of 
the Sodium layer and for the two laser projection configura- 
tions, central and lateral. The noise is clearly larger for lat- 
eral than for central projection, by a factor ~ 2. The target 
WFE in each plot is shown by the horizontal dashed-dotted 
line. Figure [14] shows the results for the bi-modal Sodium 
profile and central laser projection. In this case, the window 
used for the spot position calculations has been suitably 
optimized (see Section 13. 2p , in order to select the sharper 
component of the bi-modal spot. Comparing Figure[T3]- left 
and 1141 it may be noticed that the results are very similar, 
for a given laser projection configuration: for this reason, we 
consider in the following the Gaussian Sodium profile as the 
reference case. 

The results concerning the required number of detected 
photons to achieve the target value of WFE for different 
values of RON and in the three relevant cases mentioned 
above are gathered in Table [2] A first evident conclusion is 
that the central laser projection is more effective than the 
lateral one from the point of view of photon requirements. 
However the latter may be advantageous for other reasons, 
in particular for the background contamination due to the 
Rayleigh scattering of the up-going laser beam in the lower 
atmosphere, that is stronger in the case of central projec- 
tion. On the other hand, an additional aspect to investigate 
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RON 


Central laser 


Lateral laser 


No elongation 


le" 


380 


890 


150 


3e- 


520 


1190 


210 


5e- 


700 


1520 


270 


7e- 


910 


1910 


330 


9e- 


1100 


2200 


410 


lle- 


1270 


2400 


500 



Table 2. Number of detected photons per subaperture per frame 
required to achieve a WFE = 200nm per LGS. Three cases are 
considered: central laser projection, lateral projection and no spot 
elongation, a condition that might be achieved with a pulsed laser 
by range gating or dynamic refocusing. The results are presented 
for the Gaussian sodium profile. For the bi-modal sodium pro- 
file considered in this paper, assuming a proper windowing of the 
subaperture, the requested number of detected photons per sub- 
aperture is close to the results obtained in the Gaussian profile 
case. 



is the perspective elongation compared to the isoplanatic 
angle at the laser wavelength: the lateral launching is more 
sensitive to this problem than the central one. Finally it has 
to be noted that the considerable gain in the case without 
elongation is achievable by means of techniques such as dy- 
namic refocus, which is technologically challenging, or range 
gating, that wastes a large fraction of the laser power. 



6 CONCLUSIONS AND FUTURE WORK 

The effect of the LGS perspective elongation on the per- 
formance of a Shack-Hartmann WFS for a large aperture 
telescope has been analyzed in this paper. The Weighted 
Center of Gravity algorithm has been used for the spot po- 
sition measurement in the subapertures. The OPD measure- 
ment error in the subapertures has been propagated to WFE 
(RMS wavefront error) across the whole pupil, assuming a 
modal least square wavefront reconstruction; a fast method 
to perform this error propagation has been proposed and val- 
idated, improving the efficiency of the calculation for a sys- 
tem with a very large number of subapertures. The method 
proposed here is based on a single channel approach, i.e. 
it allows the computation of the propagated wavefront error 
per single LGS, starting from the first order properties of the 
WFS (number of subapertures, number of detected photons 
per subaperture, RON, focal plane sampling, subaperture 
FoV, calculation window). 

The analysis presented in this paper has been tuned for 
the specific case of a MCAO module for the E-ELT, assum- 
ing a target WFE per LGS of 200 nm. Different Sodium 
profiles have been considered: Gaussian and bi-modal pro- 
files. It has been found that, with a suitable optimization 
of the WCoG algorithm, performance is quite similar in the 
two cases. Two laser projection schemes have been consid- 
ered: central and lateral projection. In this framework, as- 
suming for instance a reasonable detector noise RON = 3 
e~ /pixel, it has been found that approximately Nph ~ 520 
detected photons per subaperture per frame are required 
for the central laser projection. For compariso n, the photo n 
return of the Keck LGS system () Wizinowich et al.ll2006h . 
properly adapted to the case of the MCAO system for the E- 
ELT considering all relevant parameters, would correspond 



to approximately Nph ~ 70 — 180; therefore the number 
of detected photons required to achieve our target WFE is 
a factor ~ 3 — 7.5 larger. For the lateral laser projection, 
the required number of photons is Nph = 1200, i.e. a fac- 
tor ~ 2.2 larger that for the central projection. The choice 
between central and lateral projection depends not only on 
the number of required photons, but also on other aspects 
(contamination due to the Rayleigh scattering of the laser 
light in the lower atmosphere, size of the isoplanatic an- 
gle compared to the perspective elongation), that are be- 
yond the scope of this paper. A notable reduction in the 
required number of detected photons {Nph — 210) might be 
obtained by canceling the perspective elongation by means 
of a dynamic refocus scheme combined with a pulsed laser 
(Baranec & Dekany 2008, Beckers et al. 2003, Georges et al. 
2003), which is tech nologically challenging, or range gating 
dMartin et al.l lioosl l . that however wastes a non negligible 
fraction of laser power. Other approaches are under investi- 
gation to achieve sim ilar results without fast moving parts 
jSchreiber et al.ll2008l ). 

A detailed analysis of the influence of the RON has 
been performed. Assuming a more conservative value RON 
= 5 e~ /pixel, the required number of detected photons in- 
creases by a factor ~ 1.4 with respect to the case RON — 3 
e~ /pixel. Another important aspect that has been analyzed 
is the number of pixels per subaperture, that determines 
the detector size. With typical Sodium profiles, assuming 
Nyquist sampling in the non elongated direction and central 
laser projection, the minimum number of pixels per subaper- 
ture is 12 X 12, so that 84 subapertures across the diameter 
(Table [T| corresponds to a detector size of 1000 x 1000 
pixels, that doubles in the case of lateral laser projection. 
This detector format, especially combined with a RON — 3 
e~ /pixel or better looks a demanding requirement. 

The analysis presented in this paper is based on the 
WCoG algorithm. However it is worth comparing this 
method to other algorithms (e.g. correlation maximization, 
quad-cell approach) with the aim of relaxing the require- 
ments on the number of photons and on the detector spec- 
ifications. Additional aspects to be considered in the eval- 
uation of the algorithms are the variation of the Sodium 
profile with time and the impact of the differential aber- 
rations between science and LGS channel, that usually de- 
pend on the Zenithal angle variation during the exposure. 
The Sodium profile variations imply continuous updating 
of the algorithm reference (the weighting function and the 
calibration curve in the case of the WCoG) ; the differential 
aberrations may translate into time dependent spot offsets 
in the subapertures, that need to be compensated or cali- 
brated. The evaluation of different algorithms, considering 
these additional is sues, is the subject of an additional work 
in preparation (fSchreiber et al.ll2009l ). 
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